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Abstract 

We present the first measurement of associated direct photon + muon production in 
hadronic collisions, from a sample of 1.8 TeV pp collisions recorded with the Collider 
Detector at Fermilab. Quantum chromodynamics (QCD) predicts that these events 
are primarily from the Compton scattering process eg -^ 07, with the final state charm 
quark producing a muon. Hence this measurement is sensitive to the charm quark 
content of the proton. The measured cross section of 29 it 9 pb is compared to a leading- 
order QCD parton shower model as well as a next-to-leading-order QCD calculation. 



Direct photon production in hadronic interactions has historically provided precise tests 
of quantum chromodynaniics (QCD). The associated production of charm quarks in direct 
photon events provides a potential probe of the charm quark content of the proton through 
the Compton scattering process eg — > 07 as shown in Fig. |^a [|^ . In a previous publication 
we tagged the charm quark with a full reconstruction of a D* meson. In this Letter we 
present the first measurement of the associated production cross section of 7 + /i in hadron- 
hadron collisions. In addition, since the previous publication next-to-leading order QCD 
(NLO QCD) calculations [|l| have become available. At NLO some new processes, which 
produce a significant increase in the number of photon+charm quark events, enter the QCD 
calculations. The largest of these new processes is gg -^ cc'^^ which is depicted in Fig. |I|b. 

The data for this analysis consists of an integrated luminosity of 13.2 pb^^ of "pp colli- 
sions collected by the Collider Detector at Fermilab (CDF) in the 1992-93 Tevatron collider 
run (Run la). The CDF detector has been described in detail elsewhere [Q. The events 
in the photon data sample discussed in this paper contain a cluster of energy in the cen- 
tral electromagnetic calorimeter |?7^| < 0.9, with no charged tracks pointing to the cluster. 
The clusters are required to have a transverse momentum Pt{= Psin{6)) between 17 and 
40 GeV and to be isolated, with less than 2 GeV of additional transverse energy in a cone of 



AR = a/A02 -|- Ar/2 = 0.7 around the cluster. Additional photon cuts were used which were 
identical to those used in the Run la CDF inclusive photon analysis [Q, after these cuts 
108K photon candidate events remain. Direct photon backgrounds from vr*^ and rj meson 
decays remain in the sample. They are subtracted on a statistical basis by using the photon 

6 



background subtraction "profile" method described in reference Q. Tlie upper cut on pho- 
ton Pt described above is necessary in order to use this technique. Using this method, the 
data sample is estimated to contain 35 ± 7% background from meson decays, leaving 70K 
direct photon events before muon selection. 

Muon candidates in this direct photon sample were selected by additionally requiring a 
match between a charged track with Pj- > 4 GeV in the central tracking chamber and a 
track in the appropriate muon system. For \riy\ < 0.6 the central muon system and central 
muon upgrade provide mostly redundant coverage and matching tracks were required in 
each system, while for 0.6 < |?7^| < 1.0 the track was required to be in the central muon 
extension. All three muon systems and their respective muon identification are discussed in 
detail in reference [Q, the identification is similar to that used in reference 0. After the track 
matching requirement there are 140±39 direct photon events with muon candidates (after the 
photon background subtraction described above) . The azimuthal angular difference between 
the direct photon and the muon candidate is shown in Fig. 0, showing that they are very 
back-to-back as expected for the 2-^2 scattering process. 

A number of quantities that define the quality of the muon candidates, such as the po- 
sition match between the track and the muon stub, were compared to a clean sample of 
muons from J/ijj decays, and were found to be in excellent agreement. A small background 
due to accidental track-stub matches was found in the central muon extension, and sub- 
tracted as discussed in reference 0. No corresponding background was found in the other 
muon systems. Note that some of the muon candidates are muons from charged pion and 
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kaon decays (decay-in-flight), while a smaller fraction are from charged hadrons that do not 
interact significantly in the material in front of the muon detectors (punch-through). 

These backgrounds are estimated with a study that includes both data and Monte Carlo 
simulation. Starting with the parent inclusive photon -|- jet data sample, where the main 
muon backgrounds come from, we measure the 4-vector of each charged particle. These 
charged particles include those from the jet as well as from the underlying event in the col- 
lision. We then pass this 4-vector into a CDF detector simulation as a charged pion or kaon 
or proton. Finally the results of the simulation are passed through the muon reconstruction, 
and the number of charged hadrons that faked a muon is recorded. Approximately 75% of 
the background is from decay- in- flight. This background estimate should be robust since 
the charged particle spectrum is measured directly from the parent data sample and the 
kinematics of decay-in-flight are well known. The backgrounds from non-interacting protons 
are negligible due to their relatively large interaction cross section, thus only the particle 
fractions of pions and kaons need to be determined to complete the muon background es- 
timate. A number of pion and kaon particle fraction measurements exist at LEP and the 
Tevatron [^, and they are quite consistent with one another. We have used the LEP mea- 
surements of 65% charged pion and 25% charged kaon as the nominal estimate since the 
kinematics of those measurements are best matched to the current analysis. We have used 
the other measurements to determine the systematic uncertainty (the largest deviation being 
60% pion and 20% kaon). Combining the simulation results with the particle fractions gives 
the muon background rate in the direct photon sample. Fig. ^ shows the muon Pt in the di- 



rect photon sample, along with the estimated muon backgrounds. One sees that the number 
of muon candidates is significantly larger than the estimated background. The estimated 
background can also be checked directly by the data by studying pions from reconstructed 
Kg decays. This estimate agrees well with the simulation result, but is statistically weaker, 
thus the simulation result will be used. 

Subtracting the estimated background leaves 72 ± 20 (statistical) direct photon + muon 
signal events. The photon-muon cross section is derived from the 72 photon-muon events 
by dividing by the luminosity, 13.2 p6~^, and the efficiencies for detecting the photon within 
|?7^| < 0.9 and the muon within \r]^\ < 1.0. These efficiencies are measured by a combination 
of simulation and data studies, which are described in the earlier references. The photon 
efficiency is ~ 33% with a small photon Pt dependence, while the muon efficiency is ~ 57% 
and depends slightly on the specific muon detector. The resulting photon-muon cross section 
is a'^"*" = 29 pb ±8 pb (statistical). 

There are four significant systematic uncertainties on direct photon + muon cross section: 
1) 10% from the muon background subtraction, which mostly comes from the uncertainties 
in the estimated pion and kaon particle fractions, 2) 9% from the photon background sub- 
traction uncertainty, 3) 5% from the uncertainty in the photon and muon cut efficiencies, 
and 4) 3.6% from the uncertainty in the CDF luminosity measurement. These added in 
quadrature are a 4 pb uncertainty, much smaller than the statistical uncertainty in the mea- 
surement. The measured photon-muon cross section including systematic uncertainties is 
therefore o"'^"*" = 29 pb ±9 pb (statistical+systematic). 
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The photon-muon cross section will be compared to two different QCD calculations of 
photon-muon production. The first QCD calculation is the PYTHIA [^ Monte Carlo which 
only has the leading-order contributions to the photon+heavy quark cross section, but has 
the full parton shower and fragmentation effects, and includes the photon+bottom process 
(about 25% of the total cross section). The second QCD calculation is the NLO QCD 
photon+charm calculation detailed in reference [Q. This has additional processes not present 
at leading-order, of which Fig. |I|b is the largest contribution. But with only the charm quark 
in the final state in the NLO QCD calculation, we need to convolute it with the effects 
of charm quark fragmentation. We will use the Peterson fragmentation model in PYTHIA 
to do this. An example of the resulting probability of observing a 4 GeV muon from this 
model is shown in Fig. ^ as a function of charm quark Pt- This curve includes the 10% 
branching ratio and the charm fragmentation effects. This calculation uses the massless 
quark approximation, which is adequate for photon-|-charm since the scale of the process is 
well above the charm mass, but cannot be applied to the photon+bottom process. 

The resulting QCD cross sections for photon+muon production are tabulated in table 
1, along with the measured value. While both models are consistent with the data within 
two standard deviations, the NLO QCD predictions are significantly larger than PYTHIA, 
and closer to the measured value. The ratio of the measurement to PYTHIA is very similar 
to that observed in our previous publication using reconstructed D*s to tag the charm 
quark 0. A number of variations of the NLO QCD calculation were investigated, such 
as renormalization scale, parton distributions, Peterson parameters. The largest variations 
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came from the Peterson parameters: a ±0.02 change in the fragmentation parameter caused 
the NLO QCD cross section to change by ±6%. 



Photon+Muon Data 


PYTHIA 
(CTEQ4L, fi = Pt/2) 


NLO QCD 

(CTEQ4M, /i = Pt/2) 


29 ± 9 pb 


13 pb 


21 pb 



Table 1: Comparison of the measured photon+muon cross section with the QCD models. 



In summary, the first measurement of direct photon plus associated muon production in 
hadronic interactions has been presented. Comparisons with NLO QCD photon+heavy quark 
production have also been presented, and are consistent with the measurement. 
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Figure 1: Two diagrams depicting the QCD picture of photon+charm quark production. In 
a) the dominant diagram is depicted, the Compton scattering process which enters at leading- 
order in QCD. Fig. b) represents the largest process that enters the QCD calculation for the 

first time at next-to-leading order. 
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Figure 2: The difference in azimuthal angle between the photon and the muon candidate is 
shown for the direct photon events. 
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Figure 3: The muon pt distribution in direct photon events is shown. The sum of all muon 
backgrounds is also shown. 
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Figure 4: The probability from PYTHIA for a charm quark to give a >4 GeV muon is shown 
as a function of the charm quark Pt- 
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